Considerations Around Structure-Based Drug Discovery for KRAS Using DOCK.

Molecular docking is a popular computational tool in drug discovery. Leveraging structural information, docking software predicts binding poses of small molecules to cavities on the surfaces of proteins. Virtual screening for ligand discovery is a useful application of docking software. In this chapter, using the enigmatic KRAS protein as an example system, we endeavor to teach the reader about best practices for performing molecular docking with UCSF DOCK. We discuss methods for virtual screening and docking molecules on KRAS. We present the following six points to optimize our docking setup for prosecuting a virtual screen: protein structure choice, pocket selection, optimization of the scoring function, modification of sampling spheres and sampling procedures, choosing an appropriate portion of chemical space to dock, and the choice of which top scoring molecules to pick for purchase.

Production of Isotopically Labeled KRAS4b.

Isotopically labelled proteins are important reagents in structural biology as well as in targeted drug development. The field continues to advance with complex multi-isotope labeling. We have combined our experience in high-level soluble KRAS4b expression with protocols for isotope incorporation, to achieve reliable and efficient approaches for several labeling strategies. Typical experiments achieve nearly 100% 15N incorporation, with yields in the range of 1.3-24.6 mg/L (median = 6.4 mg/L, n = 53). Improvements in the growth parameters in the presence of deuterium reduce the standard time of fermentation from 5 days to 3 days by modifying the medium used during the weaning process. The methods described are compatible with multi-isotope labeling and site-specific labeling.

Luminescence-Based Techniques for KRAS Thermal Stability Monitoring.

Various biochemical methods have been introduced to detect and characterize KRAS activity and interactions, from which the vast majority is based on luminescence detection in its varying forms. Among these methods, thermal stability assays, using luminophore-conjugated proteins or external environment sensing dyes, are widely used. In this chapter, we describe methods enabling KRAS stability monitoring in vitro, with an emphasis on ligand-induced stability. This chapter focuses mainly on luminescence-based techniques utilizing external dye molecules and fluorescence detection.

From covalent transition states in chemistry to noncovalent in biology: from ß- to Φ-value analysis of protein folding.

Solving the mechanism of a chemical reaction requires determining the structures of all the ground states on the pathway and the elusive transition states linking them. 2024 is the centenary of Brønsted's landmark paper that introduced the ß-value and structure-activity studies as the only experimental means to infer the structures of transition states. It involves making systematic small changes in the covalent structure of the reactants and analysing changes in activation and equilibrium-free energies. Protein engineering was introduced for an analogous procedure, Φ-value analysis, to analyse the noncovalent interactions in proteins central to biological chemistry. The methodology was developed first by analysing noncovalent interactions in transition states in enzyme catalysis. The mature procedure was then applied to study transition states in the pathway of protein folding - 'part (b) of the protein folding problem'. This review describes the development of Φ-value analysis of transition states and compares and contrasts the interpretation of ß- and Φ-values and their limitations. Φ-analysis afforded the first description of transition states in protein folding at the level of individual residues. It revealed the nucleation-condensation folding mechanism of protein domains with the transition state as an expanded, distorted native structure, containing little fully formed secondary structure but many weak tertiary interactions. A spectrum of transition states with various degrees of structural polarisation was then uncovered that spanned from nucleation-condensation to the framework mechanism of fully formed secondary structure. Φ-analysis revealed how movement of the expanded transition state on an energy landscape accommodates the transition from framework to nucleation-condensation mechanisms with a malleability of structure as a unifying feature of folding mechanisms. Such movement follows the rubric of analysis of classical covalent chemical mechanisms that began with Brønsted. Φ-values are used to benchmark computer simulation, and Φ and simulation combine to describe folding pathways at atomic resolution.

[Preparation of a macroporous anion exchange chromatographic medium by polyamine grafting and evaluation of its protein-adsorption behavior].

The macroporous anion exchange chromatographic medium (FastSep-PAA) was prepared through grafting polyallylamine (PAA) onto polyacrylate macroporous microspheres (FastSep-epoxy). The effects of the synthesis conditions, including the PAA concentration, reaction time, and reaction solution pH, on the ion exchange (IC) of the medium were investigated in detail. When the PAA concentration, reaction time, and reaction solution pH were increased, the IC of the medium increased, and optimal synthesis conditions were then selected in combination with changes of protein binding capacity. A scanning electron microscope was used to examine the surface morphology of the medium. The medium possessed high pore connectivity. Furthermore, the pore structure of the medium was preserved after the grafting of PAA onto the macroporous microspheres. This finding demonstrates that the density of the PAA ligands does not appear to have any discernible impact on the structure of the medium; that is, no difference in the structure of the medium is observed before and after the grafting of PAA onto the microspheres. The pore size and pore-size distribution of the medium before and after grafting were determined by mercury intrusion porosimetry and the nitrogen adsorption method to investigate the relationship between pore size (measured in the range of 300-1000 nm) and protein adsorption. When the pore size of the medium was increased, its protein binding capacity did not exhibit any substantial decrease. An increase in pore size may hasten the mass transfer of proteins within the medium. Among the media prepared, that with a pore size of 400 nm exhibited the highest dynamic-binding capacity (DBC: 70.3 g/L at 126 cm/h). The large specific surface area of the medium and its increased number of protein adsorption sites appeared to positively influence its DBC. When the flow rate was increased, the protein DBC decreased in media with original pore sizes of less than 700 nm. In the case of the medium with an original pore size of 1000 nm, the protein DBC was independent of the flow rate. The protein DBC decreased by 3.5% when the flow rate was increased from 126 to 628 cm/h. In addition, the protein DBC was maintained at 57.7 g/L even when the flow velocity was 628 cm/h. This finding reveals that the diffusion rate of protein molecules at this pore size is less restricted and that the prepared medium has excellent mass-transfer performance. These results confirm that the macroporous polymer anion exchange chromatographic medium developed in this study has great potential for the high-throughput separation of proteins.

Diverse crystalline protein scaffolds through metal-dependent polymorphism.

As protein crystals are increasingly finding diverse applications as scaffolds, controlled crystal polymorphism presents a facile strategy to form crystalline assemblies with controllable porosity with minimal to no protein engineering. Polymorphs of consensus tetratricopeptide repeat proteins with varying porosity were obtained through co-crystallization with metal salts, exploiting the innate metal ion geometric requirements. A single structurally exposed negative amino acid cluster was responsible for metal coordination, despite the abundance of negatively charged residues. Density functional theory calculations showed that while most of the crystals were the most thermodynamically stable assemblies, some were kinetically trapped states. Thus, crystalline porosity diversity is achieved and controlled with metal coordination, opening a new scope in the application of proteins as biocompatible protein-metal-organic frameworks (POFs). In addition, metal-dependent polymorphic crystals allow direct comparison of metal coordination preferences.

Investigating direct current potentials that affect native protein conformation during trapped ion mobility spectrometry-mass spectrometry.

Trapped ion mobility spectrometry-time-of-flight mass spectrometry (TIMS-TOFMS) has emerged as a tool to study protein conformational states. In TIMS, gas-phase ions are guided across the IM stages by applying direct current (DC) potentials (D1-6), which, however, might induce changes in protein structures through collisional activation. To define conditions for native protein analysis, we evaluated the influence of these DC potentials using the metalloenzyme bovine carbonic anhydrase (BCA) as primary test compound. The variation of DC potentials did not change BCA-ion charge and heme content but affected (relative) charge-state intensities and adduct retention. Constructed extracted-ion mobilograms and corresponding collisional cross-section (CCS) profiles gave useful insights in (alterations of) protein conformational state. For BCA, the D3 and D6 potential (which are applied between the deflection transfer and funnel 1 [F1] and the accumulation exit and the start of the ramp, respectively) had most profound effects, showing multimodal CCS distributions at higher potentials indicating gradual unfolding. The other DC potentials only marginally altered the CCS profiles of BCA. To allow for more general conclusions, five additional proteins of diverse molecular weight and conformational stability were analyzed, and for the main protein charge states, CCS profiles were constructed. Principal component analysis (PCA) of the obtained data showed that D1 and D3 exhibit the highest degree of correlation with the ratio of folded and unfolded protein (F/U) as extracted from the mobilograms obtained per set D potential. The correlation of D6 with F/U and protein charge were similar, and D2, D4, and D5 showed an inverse correlation with F/U but were correlated with protein charge. Although DC boundary values for induced conformational changes appeared protein dependent, a set of DC values could be determined, which assured native analysis of most proteins.

Orientational Pathways during Protein Translocation through Polymer-Modified Nanopores.

Protein translocation through nanopores holds significant promise for applications in biotechnology, biomolecular analysis, and medicine. However, the interpretation of signals generated by the translocation of the protein remains challenging. In this way, it is crucial to gain a comprehensive understanding on how macromolecules translocate through a nanopore and to identify what are the critical parameters that govern the process. In this study, we investigate the interplay between protein charge regulation, orientation, and nanopore surface modifications using a theoretical framework that allows us to explicitly take into account the acid-base reactions of the titrable amino acids in the proteins and in the polyelectrolytes grafted to the nanopore surface. Our goal is to thoroughly characterize the translocation process of different proteins (GFP, ß-lactoglobulin, lysozyme, and RNase) through nanopores modified with weak polyacids. Our calculations show that the charge regulation mechanism exerts a profound effect on the translocation process. The pH-dependent interactions between proteins and charged polymers within the nanopore lead to diverse free energy landscapes with barriers, wells, and flat regions dictating translocation efficiency. Comparison of different proteins allows us to identify the significance of protein isoelectric point, size, and morphology in the translocation behavior. Taking advantage of these insights, we propose pH-responsive nanopores that can load proteins at one pH and release them at another, offering opportunities for controlled protein delivery, separation, and sensing applications.

Differential Interactions of the Proteasome Inhibitor PI31 with Constitutive and Immuno-20S Proteasomes.

PI31 (Proteasome Inhibitor of 31,000 Da) is a 20S proteasome binding protein originally identified as an in vitro inhibitor of 20S proteasome proteolytic activity. Recently reported cryo-electron microscopy structures of 20S-PI31 complexes have revealed that the natively disordered proline-rich C-terminus of PI31 enters the central chamber in the interior of the 20S proteasome and interacts directly with the proteasome's multiple catalytic threonine residues in a manner predicted to inhibit their enzymatic function while evading its own proteolysis. Higher eukaryotes express an alternative form of the 20S proteasome (termed "immuno-proteasome") that features genetically and functionally distinct catalytic subunits. The effect of PI31 on immuno-proteasome function is unknown. We examine the relative inhibitory effects of PI31 on purified constitutive (20Sc) and immuno-(20Si) 20S proteasomes in vitro and show that PI31 inhibits 20Si hydrolytic activity to a significantly lesser degree than that of 20Sc. Unlike 20Sc, 20Si hydrolyzes the carboxyl-terminus of PI31 and this effect contributes to the reduced inhibitory activity of PI31 toward 20Si. Conversely, loss of 20Sc inhibition by PI31 point mutants leads to PI31 degradation by 20Sc. These results demonstrate unexpected differential interactions of PI31 with 20Sc and 20Si and document their functional consequences.

A new paradigm for applying deep learning to protein-ligand interaction prediction.

Protein-ligand interaction prediction presents a significant challenge in drug design. Numerous machine learning and deep learning (DL) models have been developed to accurately identify docking poses of ligands and active compounds against specific targets. However, current models often suffer from inadequate accuracy or lack practical physical significance in their scoring systems. In this research paper, we introduce IGModel, a novel approach that utilizes the geometric information of protein-ligand complexes as input for predicting the root mean square deviation of docking poses and the binding strength (pKd, the negative value of the logarithm of binding affinity) within the same prediction framework. This ensures that the output scores carry intuitive meaning. We extensively evaluate the performance of IGModel on various docking power test sets, including the CASF-2016 benchmark, PDBbind-CrossDocked-Core and DISCO set, consistently achieving state-of-the-art accuracies. Furthermore, we assess IGModel's generalizability and robustness by evaluating it on unbiased test sets and sets containing target structures generated by AlphaFold2. The exceptional performance of IGModel on these sets demonstrates its efficacy. Additionally, we visualize the latent space of protein-ligand interactions encoded by IGModel and conduct interpretability analysis, providing valuable insights. This study presents a novel framework for DL-based prediction of protein-ligand interactions, contributing to the advancement of this field. The IGModel is available at GitHub repository https://github.com/zchwang/IGModel.

SPDesign: protein sequence designer based on structural sequence profile using ultrafast shape recognition.

Protein sequence design can provide valuable insights into biopharmaceuticals and disease treatments. Currently, most protein sequence design methods based on deep learning focus on network architecture optimization, while ignoring protein-specific physicochemical features. Inspired by the successful application of structure templates and pre-trained models in the protein structure prediction, we explored whether the representation of structural sequence profile can be used for protein sequence design. In this work, we propose SPDesign, a method for protein sequence design based on structural sequence profile using ultrafast shape recognition. Given an input backbone structure, SPDesign utilizes ultrafast shape recognition vectors to accelerate the search for similar protein structures in our in-house PAcluster80 structure database and then extracts the sequence profile through structure alignment. Combined with structural pre-trained knowledge and geometric features, they are further fed into an enhanced graph neural network for sequence prediction. The results show that SPDesign significantly outperforms the state-of-the-art methods, such as ProteinMPNN, Pifold and LM-Design, leading to 21.89%, 15.54% and 11.4% accuracy gains in sequence recovery rate on CATH 4.2 benchmark, respectively. Encouraging results also have been achieved on orphan and de novo (designed) benchmarks with few homologous sequences. Furthermore, analysis conducted by the PDBench tool suggests that SPDesign performs well in subdivided structures. More interestingly, we found that SPDesign can well reconstruct the sequences of some proteins that have similar structures but different sequences. Finally, the structural modeling verification experiment indicates that the sequences designed by SPDesign can fold into the native structures more accurately.

PeSTo-Carbs: Geometric Deep Learning for Prediction of Protein-Carbohydrate Binding Interfaces.

The Protein Structure Transformer (PeSTo), a geometric transformer, has exhibited exceptional performance in predicting protein-protein binding interfaces and distinguishing interfaces with nucleic acids, lipids, small molecules, and ions. In this study, we introduce PeSTo-Carbs, an extension of PeSTo specifically engineered to predict protein-carbohydrate binding interfaces. We evaluate the performance of this approach using independent test sets and compare them with those of previous methods. Furthermore, we highlight the model's capability to specialize in predicting interfaces involving cyclodextrins, a biologically and pharmaceutically significant class of carbohydrates. Our method consistently achieves remarkable accuracy despite the scarcity of available structural data for cyclodextrins.

Structure-Based Protein Assembly Simulations Including Various Binding Sites and Conformations.

Many biological functions are mediated by large complexes formed by multiple proteins and other cellular macromolecules. Recent progress in experimental structure determination, as well as in integrative modeling and protein structure prediction using deep learning approaches, has resulted in a rapid increase in the number of solved multiprotein assemblies. However, the assembly process of large complexes from their components is much less well-studied. We introduce a rapid computational structure-based (SB) model, GoCa, that allows to follow the assembly process of large multiprotein complexes based on a known native structure. Beyond existing SB Go̅-type models, it distinguishes between intra- and intersubunit interactions, allowing us to include coupled folding and binding. It accounts automatically for the permutation of identical subunits in a complex and allows the definition of multiple minima (native) structures in the case of proteins that undergo global transitions during assembly. The model is successfully tested on several multiprotein complexes. The source code of the GoCa program including a tutorial is publicly available on Github: https://github.com/ZachariasLab/GoCa. We also provide a web source that allows users to quickly generate the necessary input files for a GoCa simulation: https://goca.t38webservices.nat.tum.de.

SadNet: a novel multimodal fusion network for protein-ligand binding affinity prediction.

Protein-ligand binding affinity prediction plays an important role in the field of drug discovery. Existing deep learning-based approaches have significantly improved the efficiency of protein-ligand binding affinity prediction through their excellent inductive bias capability. However, these methods only focus on fragmented three-dimensional data, which truncates the integrity of pocket data, leading to the neglect of potential long-range interactions. In this paper, we propose a dual-stream framework, with amino acid sequence assisting the atomic data fusion for graph neural network (termed SadNet), to fuse both 3D atomic data and sequence data for more accurate prediction results. In detail, SadNet consists of a pocket module and a sequence module. The sequence module expands the "receptive field" of the pocket module through a mid-term virtual node fusion. To better integrate sequence-level information from the sequence module and 3D structural information from the pocket module, we incorporate structural information for each amino acid within the sequence module. Besides, to better understand the intrinsic relationship between sequences and 3D atomic information, our SadNet utilizes information stacking from both the early stage and later stage. Experimental results on publicly available benchmark datasets demonstrate the superiority of the proposed dual-stream approach over the state-of-the-art alternatives. The code of this work is available online at https://github.com/wardhong/SadNet.

Single-Molecule Diffusivity Quantification Unveils Ubiquitous Net Charge-Driven Protein-Protein Interaction.

Recent microscopy and nuclear magnetic resonance (NMR) studies have noticed substantial suppression of intracellular diffusion for positively charged proteins, suggesting an overlooked role of electrostatic attraction in nonspecific protein interactions in a predominantly negatively charged intracellular environment. Utilizing single-molecule detection and statistics, here, we quantify in aqueous solutions how protein diffusion, in the limit of low diffuser concentration to avoid aggregate/coacervate formation, is modulated by differently charged interactor proteins over wide concentration ranges. We thus report substantially suppressed diffusion when oppositely charged interactors are added at parts per million levels, yet unvaried diffusivities when same-charge interactors are added beyond 1%. The electrostatic attraction-driven suppression of diffusion is sensitive to the protein net charge states, as probed by varying the solution pH and ionic strength or chemically modifying the proteins and is robust across different diffuser-interactor pairs. By converting the measured diffusivities to diffuser diameters, we further show that in the limit of excess interactors, a positively charged diffuser molecule effectively drags along just one monolayer of negatively charged interactors, where further interactions stop. We thus unveil ubiquitous, net charge-driven protein-protein interactions and shed new light on the mechanism of charge-based diffusion suppression in living cells.

G Protein Activation Occurs via a Largely Universal Mechanism.

Understanding how signaling proteins like G proteins are allosterically activated is a long-standing challenge with significant biological and medical implications. Because it is difficult to directly observe such dynamic processes, much of our understanding is based on inferences from a limited number of static snapshots of relevant protein structures, mutagenesis data, and patterns of sequence conservation. Here, we use computer simulations to directly interrogate allosteric coupling in six G protein α-subunit isoforms covering all four G protein families. To analyze this data, we introduce automated methods for inferring allosteric networks from simulation data and assessing how allostery is conserved or diverged among related protein isoforms. We find that the allosteric networks in these six G protein α subunits are largely conserved and consist of two pathways, which we call pathway-I and pathway-II. This analysis predicts that pathway-I is generally dominant over pathway-II, which we experimentally corroborate by showing that mutations to pathway-I perturb nucleotide exchange more than mutations to pathway-II. In the future, insights into unique elements of each G protein family could inform the design of isoform-specific drugs. More broadly, our tools should also be useful for studying allostery in other proteins and assessing the extent to which this allostery is conserved in related proteins.

Assessing Weak Anion Binding to Small Peptides.

Since Hofmeister's seminal studies in the late 19th century, it has been known that salts and buffers can drastically affect the properties of peptides and proteins. These Hofmeister effects can be conceived of in terms of three distinct phenomena/mechanisms: water-salt interactions that indirectly induce the salting-out of a protein by water sequestration by the salt, and direct salt-protein interactions that can either salt-in or salt-out the protein. Unfortunately, direct salt-protein interactions responsible for Hofmeister effects are weak and difficult to quantify. As such, they are frequently construed of as being nonspecific. Nevertheless, there has been considerable effort to better specify these interactions. Here, we use pentapeptides to demonstrate the utility of the H-dimension of nuclear magnetic resonance (NMR) spectroscopy to assess anion binding using N-H signal shifts. We qualify binding using these, demonstrating the upfield shifts induced by anion association and revealing how they are much larger than the corresponding downfield shifts induced by magnetic susceptibility and other ionic strength change effects. We also qualify binding in terms of how the pattern of signal shifts changes with point mutations. In general, we find that the observed upfield shifts are small compared with those induced by anion binding to amide-based hosts, and MD simulations suggest that this is so. Thus, charge-diffuse anions associate mostly with the nonpolar regions of the peptide rather than directly interacting with the amide N-H groups. These findings reveal the utility of 1H NMR spectroscopy for qualifying affinity to peptides─even when affinity constants are very low─and serve as a benchmark for using NMR spectroscopy to study anion binding to more complex systems.

DrugMGR: a deep bioactive molecule binding method to identify compounds targeting proteins.

MOTIVATION: Understanding the intermolecular interactions of ligand-target pairs is key to guiding the optimization of drug research on cancers, which can greatly mitigate overburden workloads for wet labs. Several improved computational methods have been introduced and exhibit promising performance for these identification tasks, but some pitfalls restrict their practical applications: (i) first, existing methods do not sufficiently consider how multigranular molecule representations influence interaction patterns between proteins and compounds; and (ii) second, existing methods seldom explicitly model the binding sites when an interaction occurs to enable better prediction and interpretation, which may lead to unexpected obstacles to biological researchers. RESULTS: To address these issues, we here present DrugMGR, a deep multigranular drug representation model capable of predicting binding affinities and regions for each ligand-target pair. We conduct consistent experiments on three benchmark datasets using existing methods and introduce a new specific dataset to better validate the prediction of binding sites. For practical application, target-specific compound identification tasks are also carried out to validate the capability of real-world compound screen. Moreover, the visualization of some practical interaction scenarios provides interpretable insights from the results of the predictions. The proposed DrugMGR achieves excellent overall performance in these datasets, exhibiting its advantages and merits against state-of-the-art methods. Thus, the downstream task of DrugMGR can be fine-tuned for identifying the potential compounds that target proteins for clinical treatment. AVAILABILITY AND IMPLEMENTATION: https://github.com/lixiaokun2020/DrugMGR.

SWISH-X, an Expanded Approach to Detect Cryptic Pockets in Proteins and at Protein-Protein Interfaces.

Protein-protein interactions mediate most molecular processes in the cell, offering a significant opportunity to expand the set of known druggable targets. Unfortunately, targeting these interactions can be challenging due to their typically flat and featureless interaction surfaces, which often change as the complex forms. Such surface changes may reveal hidden (cryptic) druggable pockets. Here, we analyze a set of well-characterized protein-protein interactions harboring cryptic pockets and investigate the predictive power of current computational methods. Based on our observations, we developed a new computational strategy, SWISH-X (SWISH Expanded), which combines the established cryptic pocket identification capabilities of SWISH with the rapid temperature range exploration of OPES MultiThermal. SWISH-X is able to reliably identify cryptic pockets at protein-protein interfaces while retaining its predictive power for revealing cryptic pockets in isolated proteins, such as TEM-1 ß-lactamase.

Maturation and Conformational Switching of a De Novo Designed Phase-Separating Polypeptide.

Cellular compartments formed by biomolecular condensation are widespread features of cell biology. These organelle-like assemblies compartmentalize macromolecules dynamically within the crowded intracellular environment. However, the intermolecular interactions that produce condensed droplets may also create arrested states and potentially pathological assemblies such as fibers, aggregates, and gels through droplet maturation. Protein liquid-liquid phase separation is a metastable process, so maturation may be an intrinsic property of phase-separating proteins, where nucleation of different phases or states arises in supersaturated condensates. Here, we describe the formation of both phase-separated droplets and proteinaceous fibers driven by a de novo designed polypeptide. We characterize the formation of supramolecular fibers in vitro and in bacterial cells. We show that client proteins can be targeted to the fibers in cells using a droplet-forming construct. Finally, we explore the interplay between phase separation and fiber formation of the de novo polypeptide, showing that the droplets mature with a post-translational switch to largely ß conformations, analogous to models of pathological phase separation.